Photon localization barrier can be overcome 
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Abstract 

In contradistinction to a widespread belief that the spatial localization of photons is restricted by a power-law falloff of 
the photon energy density, I.Bialynicki-BiruIa [Phys. Rev. Lett. 80, 5247 (1998)] has proved that any stronger - up to an 
almost exponential - falloff is allowed. We are showing that for certain specifically designed cylindrical one-photon states the 
localization is even better in lateral directions. If the photon state is built from the so-called focus wave mode, the falloff in 
the waist cross-section plane turns out to be quadratically exponential (Gaussian) and such strong localization persists in the 
course of propagation. 



While quantum electrodynamics (QED) underwent an 
impressive development and reached its maturity in the 
middle of the last century, one of its basic concepts - 
the photon wave function in free space - was deprived 
of such fortune. Although the photon wave function 
in coordinate representation was introduced already in 
1930 by Landau and Peierls the concept was found 
to suffer from inherent difficulties that were not over- 
come during the century (see review 0). The common 
explanation presented in textbooks (e.g., 0,0) may be 
summed up as follows: (i) no position operator exists for 
the photon, (ii) while the position wave function may be 
localized near a space-time point, the measurable quan- 
tities like the electromagnetic field vectors, energy, and 
the photodetection probability remain spread out due 
to their non-local relation with the position wave func- 
tion. However, just before the turn of the century both 
of these widely-espoused notions were disproved 01,0 
and in the new century a fresh interest in the photon lo- 
calization problem seems to have been awakened (see, 
e-g-, 0,0,0), meeting the needs of developments in 
near-field optics, cavity QED, and quantum computing. 

L Bialynicki-Birula writes 6] that the statement "even 
when the position wave function is strongly concentrated 
near the origin, the energy wave function is spread out 
over space asymptotically like r"*"/^ " (citation from 0, 
p. 638) is incorrect and that both wave functions may be 
strongly concentrated near the origin. He demonstrates, 
on one hand, that photons can be essentially better lo- 
calized in space - with an exponential falloff of the pho- 
ton energy density and the photodetection rates. On 
the other hand, he establishes - and it is even some- 
what startling that nobody has done it earlier - that 
certain localization restrictions arise out of a mathemat- 
ical property of the positive frequency solutions which 
therefore are of a universal character and apply not only 
to photon states but hold for all particles. More specifi- 
cally, it has been proven in the Letter for the case of 
spherically imploding-exploding one-photon wavepacket 
that the Paley- Wiener theorem allows even at instants 
of maximal localization only such asymptotic decrease of 
the modulus of the wave function with the radial distance 



r that is slower than the linear exponential one, i.e., any- 
thing slower than ~ exp(— Ar), where j4 is a constant. 
The latter is what the Paley- Wiener theorem says about 
a function whose Fourier spectrum contains no negative 
frequencies. 

The purpose of the present Letter is to indicate that 
one-photon wave functions of a specific type can break 
the localization restriction and exhibit the linear expo- 
nential and even faster falloff with the distance. Yet, 
there is no contradiction either with the result of Ref. 
or with the Paley- Wiener theorem, since the wave func- 
tions are cylindrical and exhibit an exceptionally strong 
localization in two dimensions out of three. 

As an introduction, we consider briefly the simplest 
case of a one-dimensional Landau-Peierls wave function 
in order to indicate how the Paley- Wiener theorem re- 
stricts the spatial localization of a photon. Then we study 
the radial falloff for three different cylindrical wave func- 
tions, using exactly the same formalism that has been 
presented in Ref. [g- Finally, the discussion of our re- 
sults allows us to refine the analysis given in Ref. . 

Let us consider a one-photon (Iph) state that corre- 
sponds to a plane-wave pulse propagating unidirection- 
ally, say, along the axis z, being polarized along a lateral 
axis (say, the x axis) 



(1) 



\lph) = / dkz f{kz)a'^{kz)\vac) 



where a'^(fcz) is the creation operator of a photon and 
f{kz) is a properly-normalized photon wave function in 
the momentum representation. Then the inverse Fourier 
transform (but including positive frequencies only!) 



H^,t)^^ I dkz f{kz)e 

dk /(fc)e*'=(^-'=*^ 
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(3) 



represents the corresponding position space wave func- 
tion of the photon in state \lph) (see, e.g., 0, p. 636). 
The modulus squared |4>(z,t)|^ gives the photon proba- 
bility density, i.e., the degree of localization along the axis 
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z (in the given case in the directions x and y any locaHza- 
tion is absent). If /(fc) differs from zero within a wide fre- 
quency band, the probabihty |$(2,t)| may be strongly 
locahzed around a point zq moving along the axis z with 
the speed of light c. However, since $(z,t), due to the 
absence of negative frequencies in the integral of Eq.(|21l, 
is nothing but a complex analytic signal, according to 
the Paley- Wiener theorem (or criterion) the asymptotic 
decrease of |$(z,t)| with the distance r = |z — zo| has 
to be weaker than ~ exp(— Ar), where A is a constant. 
All the more excluded are any finite-support functions 
in the role of ^{z^€). To conclude the introduction, let 
us notice that if the counterpropagating (with < 0) 
Fourier components are involved in Eq.Q, the Paley- 
Wiener theorem does not apply at the instant t — 0. 

Following Ref. 0], we shall study the photon localiza- 
tion by examining the asymptotic behavior of the positive 
frequency part of the Riemann-Silberstein vector F(r, i) 
(called the energy wave function) which directly deter- 
mines the energy density of a one-photon state and is con- 
veniently expressed through a " superpotential" Z(r, r) 



F(r,i) = V X 



V X Z(r,r) 



(4) 



where t = ct. The vector field Z(r, r) is nothing but 
an analytic signal version of the Hertz potential, i.e., 
Z(r,T) = m^(r, r), where m is a constant vector that 
includes the proper normalization factor and ^'(r, r) is 
any solution of the scalar wave equation, which is taken 
in the form of the analytic signal. 

As the first example leading to a stronger localization 
that one might expect from the Paley- Wiener theorem, 
let us consider the photon field where m is directed along 
the axis z (any other orientation gives similar results) 
and 5'(r,T) is a superposition of cylindrically symmetric 
Bessel functions Jo as a wavepacket with the exponen- 
tial spectrum and a specific dispersion law for the axial 
wavenumber kz{uj) = const — fco 



4'(p,z,t) 



dk Jo (kpp) e 



- kA —i{kT — ko z) 



(5) 



where the radial coordinate p has been introduced and 

U2 ,.2\l/2 



is the lateral component of the wave 



vector of the monochromatic plane-wave constituents 
represented with the weight function e~'^^ whose width 
is A~^. The integral can be taken with the help of a 
Laplace transform table and we obtain 



exp 



Z(p,2,r) = m- 



1^0 1 



(A -I- ir) 



ikoz 



+ + irf 



(6) 



Eq.(|SJl describes a simple cylindrical pulse modulated 
harmonically in the axial direction and radially converg- 
ing (when r < 0) to the axis and thereafter (when r > 0) 



expanding from it, the intensity distribution resembling 
an infinitely long tube coaxial with the z axis and with 
a time-dependent diameter (see Fig. 6. in Ref. ^3)- 
follows from Eqs. ® and that 



Z(p CXD, Z, T 
|F(p ^ 00,Z,T 



0)1 ^p-iexp(-p/0 



(7) 



0)1^ 



+ 0(p-3)] exp(-2p/0 



(8) 



where / = |fco|~ is the characteristic length (or length 
unit). Thus, while the photon is delocalized in the axial 
direction, its energy density falloff in the lateral direc- 
tions is exactly the linear exponential one at all times the 
conditions r ^ p A are fulfilled, see Fig. 1. The time 
derivative as well as the spatial derivatives contain the 
same exponential factor, ensuring the exponential falloff 
of the the Riemann-Silberstein vector in Eq.©. Hence, 




FIG. 1: Curves of the radial dependence in a decimal logarith- 
mic scale. Curve Ais for |Z(p,0,r = 0)|; B, |Z(p,0,r = 2.501; 
C is the same as B but with 'I' taken from Eq. @; D, 
I T^Z(p, 0, r = 0) I ; E is a reference curve exp(— p/i). The 
curves A, B, and C have been normalized so that |Z(0, 0, 0)| = 
1. The values of the remaining free parameters are A = O.IZ 
and (3 = 0.8. 



a one-photon field given by Eq. © serves as the first and 
simplest example where the localization in two transver- 
sal dimensions is governed by different rules than local- 
ization in three dimensions according to Ref. 0. 

The next example is readily available via the Lorentz 
transformation of the wave function given by Eq. (jsj along 
the axis z, which gives another possible solution of the 
scalar wave equation. The result is a new independent 
solution but it can also be considered as the wave given 
by Eqs. (O and lO, which is observed in another inertial 
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reference frame |lo| : 



^'(p,2:,r) 



\ko\ ^p^ + {A-ij {f3z-T)f'^ 
X exp (ijko {z ~ Pt)) , (9) 



where the relativistic factors 7 = (1 — P"^)"^^"^ and /3 = 
w/c < 1 have been introduced, v being a free parameter ~ 
the relative speed between the frames. In the waist region 
(see Fig. 2) this wave function has the same radial falloff 
as was given by Eq. Q , see curve " C" in Fig. 1 , while the 
axial localization follows a power law. The strongly lo- 
calized waist and the whole amplitude distribution move 
rigidly and without any spread along the axis z with a 
superluminal speed c//3. Such wave with intriguing prop- 
erties, named the focused X wave (FXW) belongs 
to the so-called propagation-invariant localized solutions 
to the wave equation - a research field emerged in the 
1980-ies (see reviews and recently reached its 

first experimental results |lfil|-[2ll|. Hence, in its waist 



(a) 



(b) 




FIG. 2: The superluminal FXW given by Eq. ©. Shown 
are the dependences (a) of the modulus and (b) of the real 
part of the wavefunction on the longitudinal {z, increasing to 
the right) and a transverse (say, x) coordinates. The distance 
between the grid lines on the basal plane {x, z) is 22A, where 
A = 27r|fco|~^, fco being negative. The values of the remaining 
free parameters are A = 30A and j3 — 0.995 or 7 = 10. 



(cross-sectional) plane a one-photon field given by the 
FXW possesses the same strong localization at any time 
as the previously considered cylindrical field does in any 
transversal plane at the instant t — 0. 

By making use of the historically first representa- 
tive of localized waves - the so-called focus wave mode 
(FWM) Hl-lil (see also Ref. To^ and reviews [13, ^ 
and references therein) one readily obtains an example of 



the field that exhibits even much stronger than exponen- 
tial localization. FWM is given by the scalar function 



*(p,z,t) 



exp 



2/(a-i(z-r)) 



i {z — t) 



exp 



i{z + t) 
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, (10) 



where again I is a wavelength-type characteristic length 
and the constant a controls the axial localization length. 
Since the FXW in the limit ^ 1 becomes a FWM 10], 
Fig. 2 gives also an idea how a FWM looks like. Multiply- 
ing Ea. Hl()(l by m to build the vector Z{p,z, t) and insert- 
ing the latter into Eq.l^J we obtain that in this example 
the photon localization in the waist plane is quadratically 
exponential (Gaussian falloff): 

|Z(p ^ cx),2 = t)| - exp(-pV2?a), (H) 
\F{p ^ 00, z = t)P - exp(-pV;a) . (12) 

In Eq. (|12|l only the highest-power term with respect to p 
is shown. 

To start discussing our results let us ask first whether 
the wave functions considered are something extraor- 
dinary. The answer is: yes, they are indeed, since 
the browsing of various integral transform tables reveals 
rather few examples where both the real and imaginary 
part of a wave function and of its time derivative have 
simultaneously an exponential or stronger localization in 
conjunction with other requisite properties. Still, the list 
of proper wave functions with an extraordinary strong lo- 
calization is not poor - in addition to an optically feasible 
version [2^ of the FWM various new interesting solutions 
can be derived [i^. Yet, it could be argued that the well- 
known Gaussian beam pulse has the same quadratically 
exponential radial profile in the waist region. However, 
resorting to the family of the Gaussian beams (the Gauss- 
Laguerre and Gauss-Hermite beams, etc.) is irrelevant 
here. The reason is that all these beams are solutions of 
the wave equation only in the paraxial approximation not 
valid in the case of any significant localization of wide- 
band (pulsed) superpositions of the beams, whereas in 
fact, e.g., an exact solution corresponding to a lowest- 
order (axisymmetric) Gaussian beam has a weak power- 
law radial falloff in the waist region [27l|.[28j|. 

The next possible objection to the physical significance 
of the results obtained might arise from the infinite to- 
tal energy of the waves given by Eqs. ©, and 
However, at any spatial location the wave func- 
tion is square integrable with respect to time, thus the 
condition of the Paley- Wiener theorem has been satis- 
fied. Moreover, physically feasible finite-energy versions 
of localized waves generally exhibit even better local- 
ization properties, although not persistently. A finite- 
energy version of the FXW, called the modified focused 
X wave (MFXW 113), has the same exponential factor as 
in Eq. (|5J| , which is multiplied by a fraction that allows to 
force the axial localization to follow an arbitrarily strong 
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power-law. The latter circumstance indicates that the 
strong lateral localization of the fields considered does 
not appear somehow at the expense of their axial local- 
ization. As a matter of fact, energy- normalization of a 
wave function depends on how many photons it describes. 
It is easy to see that our derivation and results hold for 
any number state with N > 1 and also for incoherent 
mixtures of such states (which is important for experi- 
mental studies). In contrast, the localization of coherent 
states is not governed by the Paley- Wiener theorem . 

The final crucial question is, are our results not in 
contradiction with those of Ref. 0? The answer is no, 
since in the case of the cylindrical waves the radial dis- 
tance and temporal frequency are not directly Fourier- 
conjugated variables. In order to clarify this point, let 
us first take a closer look at the proof of the Paley- 
Wiener limit for three-dimensional isotropic localization. 
In Ref. Eq.(24) for Z(r,r) contains a superposition of 
spherically symmetric standing waves 



dk h{lk) 



sin kr 



— ikr 



9i- 



(13) 

where h(lk) is the spectrum and g{.) is its Fourier image. 
The sine in Ea. (|13(l results from the imploding and ex- 
ploding spherical wave constituents of the standing wave, 
like an odd one-dimensional standing wave arises from 
counterpropagating waves. While the asymptotic behav- 
ior of the function g{.) and hence of the function Z(r, r) 
for large values of the radial distance r are generally re- 
stricted by the Paley- Wiener theorem, strictly at the in- 
stant r = of maximal localization the integral is nothing 
but the sine transform for which the theorem does not 
apply. Indeed, the sine transform tables give examples 
of the resultant functions with arbitrarily abrupt falloff. 
However, it does not mean as if the localization restric- 
tion was lifted at the instant t — 0. The explanation 
is that according to Eq.Q the energy wave function in- 
volves also the time derivative of Z(r,T), but the sine 
transforms of two functions h{lk) and h{lk)k cannot si- 
multaneously possess arbitrarily abrupt falloffs. In con- 
trast, the time derivative of the wave function given by 
Eq. © or Eq. Q has the same strong exponential falloff as 
the function itself, which persists for some (not too long) 
time, see Fig.l. By comparing Eqs. jSJ and H13|) we notice 
that while in Ea. (|13|l - as well as in its one-dimensional 
equivalent - the argument of the sine function is the prod- 
uct of the distance with the Fourier variable, in Eq. ^ the 
argument of the Bessel function is the product of the ra- 
dial distance p with the radial wavenumber kp the latter 
depending on the Fourier variable through the square- 
root expression with the constant parameter fco ~ the 
lower limit of the integration. As it follows also from 
Eqs. 10 and the condition fco 7^ is crucial for ob- 
taining the exponential falloff. Hence, in the case of the 
cylindrical waves considered by us, the apparent viola- 



tion of the rules set by the Paley- Wiener theorem results 
from the specific complicated relation between the radial 
distance and the Fourier variable. 

In conclusion, we have shown that for certain cylin- 
drical iV-photon states {N = 1,2, ...) the localization in 
lateral directions breaks the limit established in Ref. 0| 
for the case of uniform spherical wave functions. These 
results hold not only for photons but for number states 
of any particles. 

The first author is thankful to Iwo Bialynicki-Birula for 
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